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ABSTRACT 

This paper measured unsteady temperature fields of 
uncoated-monolith and catalytic monolith under real 
engine operating conditions using themiocouples. A 
multi-dimensional flow mathematical model of the 
turbulence, heat and mass transfer, and chemical 
reactions in monoliths was established using a 
computational fluid dynamics (CFD) code and 
numerically solved in the whole flow field of the catalytic 
converter. The purpose of this paper is to study 
unsteady wam-up charactBristics of the monoliths and 
to investigate effects of inlet cone structure on 
temperature distribution of the catalytic converter. 
Experimental results show that the warm-up^behaVlors 
between uncoated-monolith and catalytic monolith are 
quite different. Simulation results indicate that the 
established model can qualitatively predict the wamn-up 
characteristics. Increasing the inlet cone angle can 
improve the light-ofF characteristics of the catalysts due 
to high flow velocity and high temperature in the center 
of the monoliths. 

INTRODUCTION 

New vehicle emission regulations equivalent to Euro I 
have been implemented in China since Jan. 2000 and 
Euro 11 will also be carried out in 2004 early or later. 
Emission regulations in China should be the same as 
those developed European countries near 2010 year. 
According to existed experiences in the world the 
main emission control technology for gasoline vehicles 
which can meet Euro II and Euro HI is still the three-way 
catalytic converter (TWC) plus electronic fuel injection 
(EFI) system. But a match of TWC with EFI and an 
optimum design of catalytic converter exhaust system 
should be more careful and refined. The main challenge 
is how to reduce the hydrocarbon (HC) emitted from 
vehicles during the cold start. According to vehicle 
driving cycle tests, approximately 80 percent of HC is 
emitted in cold start phase due to a low temperature of 
the exhaust gas and monolith'^. Therefore, how to 
increase HC conversion efficiency during the cold start is 
a key factor to satisfy the stringent emission regulations. 
In this paper, temperature fields of the monoliths were 
measured and simulated to study temperature 



distributions and warm-up characteristics of the catalytic 
converter. 

MEASUREMENT OF TEMPERATURE IN MONOLITH 
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Figure 1. Schematic of test rig 




Figure 2. Installation of thermocouples in converter 



A schematic of test rig and the installation of 
thermocouples in the catalytic converter are shown in 
Figure 1 and Figure 2 respectively. The test engine is 
BN492 S.I. engine with a single point injection (SPI) 
system. Three-way catalyst monolith is used as a 
catalytic monolitii. Uncoated-monolrth is made of 
cordierite (2MgO-2Al203*5Si02) , with a porosity of 300 
cpsi. The NiCr-NtSi themnocouple used In the 
experiment is 1 mm in diameter. From the rear end of 
the monoliflis, the thermocouples are inserted into the 
channels at different positions along radial and axial 
directions (also see Figure 3 below). The inlets of the 
test channels are sealed by temperature-resisting 
cement to measure the temperature of the monoliths. 
Severe! themnocouples protrude their heads 10 mm out 
of tfie channels to measure the inflow gas temperature, 
as shown in Figure 2. Exhaust gas from the engine has 
two routine lines (1 and 2 as shown in Figure 1) to enter 
the catalytic converter, one is to pass through the heat 
exchanger, and anotiier is to enter the converter directty. 
The heat exchanger installed in routine 1 is used to cool 
the hot exhaust gas to a fixed temperature for an initial 
test temperBtJre. which is greatiy lower than the light-off 
temperature of the catalyst. The control valve turns the 
exhaust gas from routine 1 to routine 2 when the 
measurement begins. 

ARRANGEMENT OF MEASURING POINTS 
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(a) Uncoated-monolltti 
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(b) Catalytic monolith 
Rgure 3. Arrangement of measuring points 

The measuring points of the themnocoupies in 
axisymmetric planes of uncoated-^nonolith and catalytic 
monolitii are arranged as shown in Figure 3(a) and (b) 
respectively. The points are mainly distributed in the 
center and near the edge of the monoliths. Along axial 



direction, the thermocouples are located at three cross- 
sections with X coordinates of 0.1 5L, 0.5L and 0.75L 
Along radial direction, they are located at the positions 
with r coordinates of 0.1 R, 0.5R, 0.9R and 1R for 
uncoated-monolith, and 0.1R, 0.6R, 0.8R and 0.9R for 
catalytic monolith. In order to obtain the gas 
temperatures into the monolitiis, three points with r 
coordinates of 0.1 R, O.SR and 0.8R are chosen for 
uncoated-monolith, and two points with r coordinates of 
0.1 R and O.SR for catalytic monolith. 

TEST RESULTS AND ANALYSIS 

Warm-up characteristics of uncoated-monolith 

Researches on warmnjp characteristics of uncoated- 
monolitti can help us to understand tiie effect of the. 
monolith itself on light-off performance of the catalyst 
The catalytic converter used in ttiis experiment has an 
inlet cone angle of 40** . The engine running conditions 
are as follows: engine speed is 2000 rAnin; A/F (air/fuel 
ratio) is controlled near stoichiometric ratio; mass flow in 
the catalytic converter is 33.3 g/s. Befbre measuring, the 
exhaust gas passes through routine 1 pipe (see also 
Figure 1) and is cooled by the heat exchanger and then 
goes into the catalytic converter. In this case, the inlet 
gas temperature is controlled at about 150 ""C. When the 
engine runs steadily, turn the valve rapidly to routine 2 to 
let the exhaust gas enter the converter directly. At the 
same time, record the temperature rising process at 
different time using data sampling and analyzing system. 
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Figure 4. Temperature contours in uncoated-monolith at 
different time 



Figure 4{ah'(f) show temperature contour changes In 
uncoated-monolith. These contours illustrate the 
distribution and development of the temperature field in 
the monolith at x=0.15L~075L along axial direction and 
r=0.1R^1R along radial direction. The front center of the 
monolith has a highest temperature due to a high flow 
velocity in the center of the monolith ^. As time goes on, 
the high temperature field expands to the surrounding 
area and also moves backwards with the heat flux 
transfer along axial and radial directions. And the 
temperature gradient decreases gradually and the 
temperature tends to be more and more unifomn. Figure 
4(e) and (Q are temperature contours at 60 second and 
90 second respectively. From the two f^ures we can see 
that the temperature difference between the two 
moments is not obvious, which indicates that the 
temperature in uncoated-monolith tends to be stable 
after 60 second. But the temperature in the edge of the 
monolith is still lower than that in the center because of 
the heat exchange with circumference near the edge. 

Warm-up characteristics of catalytic monolith 

For catalytic monolith, not only is there a convective heat 
transfer between the monolith and the exhaust gas, but 
also a large amount of heat released from the chemical 
reactions on the surfece of the monolith. Consequently, 
the unsteady temperature field of catalytic monolith is 
quite different from uncoated-monolith. Here, all the 
engine test conditions and measuring process are the 
same as mentioned above in uncoated-monolith. But 
catalytic monolith has 15 measuring points more than 
uncoated-monolith. 
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Figure 5. Warm-up behaviors at point 1 and point 16 

Figure 5 shows the wamn-up behaviors of measuring 
point 1 and point 16, which gives the temperature of the 
front center of the monolith and the temperature of the 
exhaust gas into the monolith respectively as shown in 
Figure 3. The f^ure indicates that the temperature of the 
monolith (point 1) is lower than that of the exhaust gas 
(point 16) at the beginning of the warm-up. But after 
about 12 second, the temperature of point 1 begins to 
exceed the temperature of point 16, which means that 
the catalyst is activated and the released heat is now 
heating the monolith . 




(d) 20 second 




(f) 90 second 



Figure 6. Temperature contours in catalytic monofith at 
different time 

Figure 6(a)-'(f) are temperature contours in catalytic 
rncanolrthi atcdlffcsmit Unm. CduiiijuuK&ufa wrthi fFj^^is 4i, it isB 
found that the developnrieiit trend of unsteady 
temperature fields of cataiyUc monolith is similar to 
uncoated-monolith at the beginning (before the catalyst 
is activated), that is to say, the front center of the 
monolith has the highest temperature and then the high 
temperature field expands and moves backwards as 
time goes on. This is because the monolith is mainly 
heated by the exhaust gas during this period, and its 
temperature increases slowly. After the catalyst is 
activated, the heat from chemical reactions resute in a 
sudden rise of the temperature in the front part As a 
result, the developmerrt of temperature fieU in catalytic 
monolith is quite different from that in uncoated-monolith. 
From Figure 6(e) and (f) it can be seen that at 60 second 
the temperature in the rear center of the nrK)nolSh is 
close to the front part; but at 90 second the temperature 
in the rear center is higher than that in the front part The 
main reason is that after the catalyst is activated, the 
monolith temperature in the front part is greatly higher 
than the exhaust gas. At this time the exhaust gas 
actuail)^ Qaala>thft; monolithi. boj noiionl^ isbthe^Qhemicali 
reaction heat biOMght away by the exhaust gas, but also 
Ihe girm^ trale ^ ithe tompess^^ liia ililiie IsDist ipari ds 
ioweied by the convective heat transfer. How^iar, 
temperatures of the middle and the rear part are a little 
low at the beginning, then the hot exhaust gas and 
released chemical reaction heat from the upstream 



accelerate chemical reactions downstream. Therefore, a 
high temperature field appears in the middle and rear 
part of the monolith. 

The contours also indicate that the temperature gradient 
on the whole monolith is very small at 5 second, the 
difference between the highest temperature and the 
lowest one is about 70 X:. After the catalyst is activated, 
the difference grows qukMy due to the beat fnm the 
chemical reacSons m the montM}, for exanopfe, the 
difference at 10 second is about UOt: and at 15 second 
it is up to about 185X:; 20 second tater, tt\e difference 
begins to decrease and tends to be constant about 70X: 
as time goes on. 

NUMERICAL SIMULATION OF THE CATALYTIC 
CONVERTERS 



The gas flow in parallel channels of the monolith can be 
regarded as a fully developed laminar flow whose 
governing equation is given by Hagen-Poiseuill equation 
as follows: 

Bp -32upU 

where JJ is the gas velocity, a is the monolith porosity, 
^ 6 the hycfraulfc diametdr of monofitfi channells. The 
gas viscosfty o which rs a funcGbn of gas temperature 
is given by : 

V = (6.542X lO-"7;) +(6. 108x 10-*r,)-0.89x 10"' (2) 

Heat and mass transfer eouations 
Gas enthalpy conservation equation 



In the past people designed the catalytic converter 
mainly relying on their experiences, which not only 
wasted a fcarge amount of time but also cost much of 
money and manpower With the rapid development of 
the computer technology and CFD in recent years, 
people can \Mze rNimericsl simijtaSon nTelhod to s#udy 
the catalytic converter perfomnance ^"^l This paper uses a 
CFD code STAR-CD to simulate the transient 
temperature fields of the catalytic converter and 
investigates the effect of inlet cone structure on 
temperature distribution and wamrvup characteristics in 
monoliths. 

MODELS OF FLOW. HEAT AND MASS TRANSFER, 
CHEMICAL REACTIONS IN MONOLITHS 

A two-dimensional computational region of the axis- 
symmetric catalytic converter (also see F^ure 2) » 
divided into a free stream region and a monolith region 
as shown in Rgure 7. The monolith region can also be 
subdivided into a fluid region and a corresponding solid 
region. This paper considers the heat and mass transfer, 
chemical reactions between the gas and the monolith 
surfece, and also the heat conduction in the monolith. 
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Figure 7. Computational region of the catalytic converter 



Flowaovemino equation 

For the fluid dynamics that occur in the free stream 
region the Reynolds averaged Navier-Stockes equations 
are solved, where the standard k-s turbulence model 
is used to close the equations. In the monolith an 
equivalent continuum approach which views the 
monolith as a porous medium through v^ich the 
unidirectional gas continually passes, is utilized to 
establish the models. 



where (s the gas density; h is the gas heat transfer 
coefficient; is the ratio of monolrtti surface to monolith 
vofume; 7^ es ttre mmM\ temperafure and ffie gas 
temperature Is obtained by: 

T^^KIC^ (4) 
where \ is the gas enthalpy and is the gas specific 
heat capacity. 

Chemical species consen/ation equation 

^+£(P,t/c,) = p./:,a,(c,-cJ (5) 

where is the mass fraction of species i in the gas 
and is the mass fraction of species i on tiie surface 
of the monolrth; is the mass transfer coefficient of 
species i . The concentrations of the species on the 
catalyst surface are governed by expressions of the form: 

where is the molecular mass of species / ; is ttie 

chemical reaction rate of species 7 . The tenn on the left 
side of equation (6) represents the diffusion process of 
species to and from the monofrth surface. The temi on 
the right side indicates the removal of species by 
reaction. 

Monolith heat conduction equation 

The equation governing the monolith temperature 
behavior is essentially that of heat conduction in the 
monolitii. Considering the orthotropic nature of the heat 
conduction, the governing equation takes the form: 

a-aKP,-^-^ja-a)-^+G(^+-^)]=5.+5, (7) 

where x is the gas flow direction; y and z are anotfier 
two coordinate directions perpendicular to x direction; 



is the monolith density; is the monolith specific 
heat; A:, Is the monolith themnal conductivity; G is the 
monolrth anisotropic conductfvrty factor in y and z 
directions. The two source items and jj, are given by: 
s,^haXT,'T;) (8) 

where 5, is the heat transfer between the monolith and 
the gas; 5, is the heat released from chemical reactions; 
tsH, is the reaction heat of species i . 

Qh^miCf?<l taction ^qM^feng 

The chemical reactions of the catalysts on the monolith 
surface are very complex. This paper only considers the 
oxidation reactions of CO and C3H6 occunred on the 
surface as a preliminary invest^ation. Propylene (CaHe) 
is assumed to be a representative of liast oxidizing 
hydrocart>ons". The simplified chemical reactions can be 
expressed by: 

CO+0.5O, -►CO, (10) 
C3H.+4.5O, ->3CO,+3H,0 (11) 
and the reaction rates are given by: 



(12) 



(13) 



where 

k, = 6.699X 10' exp(-12556/7;) ( 1 5) 

it, =1.392x10" e3q)(-14556/rj (16) 

a:, =65.5xexp(961/7;) (17) 

a:, =2.08x10' exp(361/7;) (18) 

K, =3.98xexp01611/rj (19) 

K, = 4.79x 10* exp(-3733/r,) (20) 

From the stoichiometry the reaction rate for oxygen must 

satisfy the following relation: 

i?^ =0.5/2^+45;?^, (21) 

SOLUTION CONDITIONS 

The PISO algorithm is used to solve the above 
governing equations for transient flows. Subroutines for 
heat and mass transfer and chemical reactions on the 
monolith surface are programmed by authors and solved 
coupling with the flow governing equations in free stream 
region (see also Figure 7). The detailed solution process 
can be refenred to the literature [7]. 

Boundary conditions 

Inlet boundary 

The inlet velocity at the entrance of the catalytic 
converter is assumed to be unrfbrm and flows along the 



symmetrical axis of the converter. Here the inlet velocity 
is 21 .2 m/s calculated from the above measurement 
of the mass flow, and the density is assumed to be 
1.0 kgAn^. The inlet mass fractions of the spades are 
supposed as follows: c^^ =0.0048 kg/kg. c^^ =0.022 
kg*g. Cc^^ =0 .00065 kg/kg. The turbulent Wnetic energy 

K is assumed to be 0.02 and its dissipation rete 
f =0.01 mV. The inlet temperature changes as shown 
in Figure 8. 
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Figure 8. Inlet temperature vs. time 
Ouflet boundary 

The outlet boundary condition is treated as a fully 
devebped flow, namely the grediente of all variables 
(except pressure) are zero. The outiet mass flux must 
satisfy the mass consen^ation. 

Wall boundary 

Wall boundary conditions are no-slip velocity boundary 
and adiabatic temperature boundary. The inlet and outlet 
ends of the monolith are also treated as an adiabatic 
temperature boundary. 

Initial conditions 

In the whole region of the catelytic converter, the initial 
velocity and concentrations of O2, CO, C3H6 are set to 
be zero and tiie temperature is assumed to be 400 K. 

SIMULATION RESULTS AND ANALYSIS 

Warm-UD characteristics of the catalytic converter 

Figure 9(a)~(h) show the temperature contours in the 
catelytic converter with a 40^ cone angle at different 
time. The upper part of the figures illustrates the 
distribution of the monolith temperature Ts and tiie down 
part is the distribution of the gas temperature Tg. From 
the figures we can see that tiie inlet temperature of the 
monoiitti is only 500 K at 10 second and the catalyst is 
not activated. At this time the temperature behind the 
monoiitti maintains the initial temperature of 400 K. At 20 
second the inlet temperature of tiie monolith rises to 600 
K and tiie catelyst is lighted off. Then the temperature in 
the front part of the monolith increases quickly and the 
gas temperature behind the monolith also b^ins to rise 



gradually. At 30 second the inlet temperature goes up to 
700 K and the highest temperature in the front part of the 
monolith reaches 734 K due to the heat released from 
chemical reactions, which is higher than the inlet 
temperature at the moment From 30 second to 40 
second, the rising rate of the temperature in the monolith 
closes to the h^hest about 5.7 K/s. After 40 second the 
temperature changing trend becomes slow. At 60 
second the gas temperature behind the monolith is close 
to the inlet temperature. At 80 second, the temperature 
behind the monolith greatly exceeds the Inlet 
temperature and the highest tempeiature reaches 755 K. 
At the same time the temperature in the whole region of 
the catalytic converter gradually tends to be stable. At 
100 second all of the temperatures are nearly unifomi 
and stable. The highest temperature is up to 81 3 K. 
From the above analysis of the warm-up behavior in the 
catalytic converter, it can be concluded that the monolith 
explores a series of heating processes by inlet gas flow 
and chemical reactions. The temperature difference in 
the monolith changes from a small value at the 
beginning to a large one and again to a small one at the 
end. The computed wann-up characteristics of the 
monolith have a qualitative agreement with the test 
results above. This means that the established model in 
this paper can essentially display the features of the flow, 
heat and mass transfer and chemical reactions occurred 
In catalytic converters. 
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(d) 30 second 



(h) 100 second 



Figures. Temperature contours of 40** catalytic converter at 
different time (K) 



Effect of inlet cone structure on wamn-uo characteristics 

In order to investigate the effect of the inlet cone 
structure on the temperature distribution, the flow 
including heat and mass transfer and chemical reactions 
in the catalytic converter with a 1 1Q"" cone angle is also 
simulated. Compared with the 40'' catalytic converter, 
the 120*" catalytic converter has tiie same computational 
conditions except for a different cone angle. 
Figure 10{a)~(h) show the temperature contours in the 
catalytic converter with a MO"* cone ar^le at different 
time. From the comparison between Figure 9 and Figure 
10 it is found that there are some differences existed in 
the temperature distribution and wamvijp behavior 
between two converters although they have a similar 
temperature change trend. For the 120 catalytic 
converter the catalyst is not lighted off at 10 second. In 
this case two catalytic converters have no obvious 
differences for the temperature distribution along axial 
direction, but the 40 catalytic converter has a more 
unifomi temperature distribution along radial direction. 
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Figure 10. Temperature contours of 120* catalytic converter 
at different time (K) 

Figure 11 shows the comparison of the temperature 
distribution along radial direction in the middle sections 
of the monoliths at 10 second. From the figure It can be 
seen that for the 120 catalytic converter the 
temperature diff^erence between the center and the edge 
in the monolith is about 30 but for the 40'' catalytic 
converter it is less than 10 K. This temperature 
distribution is consistent with the velocity distribution in 
the same location at the same time, as shown in Figure 
12, which means that tiie gas velocity has a significant 
influence on the temperature distribution in catalytic 
converters. 
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Figure 1 1 . Temperature distribution in middle cross section 
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Figure 12. Axial velocity distribution in middle cross section 

Figure 13 and Figure 14 show the warm-up behaviors at 
the center and at the edge in the middle section of the 
monolittis respectively. From tiiese figures we can 
clearly see that the temperature of 120^ converter 
exceeds the temperature of 40* converter as time goes 
on, and the largest temperature difference is more than 
50 K at about 40 second. However, as the temperature 
tends to be stable, the difference becomes smaller and 



smaller. In addition, the difference is larger at the center 
than at the edge. This Is because the flow velocity at the 
center of 120" converter is higher than that of the 40"* 
converter (also see Figure 12) and the Inlet temperature 
can quickly conduct into the monolith, which results in a 
rapid light-oflf of the catalysts. Moreover, the heat 
released from chemical reactions accelerates the warm- 
up of the monolith. 




Figure 13. Warm-up characteristics at the center 



goes on. the high temperature region expands and 
moves backwards. 

For catalytic monolith, when the catalyst is not 
activated, the temperature distribution and warm-up 
characteristics of the monolith are similar to uncoated 
monolith. However, after the catalyst is activated, the 
high temperature appears In the rear center part of 
the monolith. 

A multi-dimensional flow including heat and mass 
transfer, chemical reactions is numerically simulated. 
Although the model contains many simplifications and 
assumptions it has been demonstrated that it can 
predict the wamrvup behavior of tiie catalytic 
converter and the effect of the inlet cone structure on 
the warm-up behavior. 

Simulation results show that the 120 catalytic 
converter has a t)etter l|ght-off perfbnnance than the 
40*" catalytic converter, especially in the center of the 
monolith. 

Further study includes quantitative validation of ttie 
model, the modeling of more tiian two reacting 
species, external heat k)ss and the theimal inertia 
effects of the catalyst matting and can. 




Figure 14. Warm-up chafactefistics a&the edge 



CONCLUSION 

From the experiment and numerical simulation of the 
temperature fields of the catalytic conve/teis, the 
conclusions can be drawn as follows: 

1 . Heat ccxwec^e Irans^ between the ei^ausl ^as 
and the monolith, heat conduction in the monolith, 
heat from chemical reactiorts, and flow velocity 
distrit)ution are all main factors affecting the 
temperature distilbution and warm-up characteristics. 

2. For uncoated-monolith, high temperature always 
appears in the front center of the monolith. As time 
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